Astronomy Sz Astrophysics manuscript no. enrich8 


February 2, 2008 


(DOI: will be inserted by hand later) 





Metal Enrichment Processes in the Intra-Cluster Medium 

S. Schindler 1 , W. Kapferer 1 , W. Domainko 1 , M. Mair 1 , E. van Kampen 1 , T. Kronberger 1 , S. 
Kimeswenger 1 , M. Ruffert 2 , O. Mangete 2 , D. Breitschwerdt 3 



1 Institut fur Astrophysik, Leopold- Franzens Universitat Innsbruck, Technikerstrafie 25, A-6020 Innsbruck, Austria 
http:/ /astro. uibk.ac.at/astroneu/hydroskiteam/index. htm 



in 
o 
o 

< 



> 
00 

O 

O 

in 
o 

S3 

Oh 
6 



X 



z School of Mathematics, University of Edinburgh, Edinburgh EH9 3JZ, Scotland, UK 
3 Institut fur Astronomie, Universitat Wien, Tiirkenschanzstr. 17, 1180 Vienna, Austria 

Received / Accepted 

Abstract. We present numerical simulations of galaxy clusters which include interaction processes between the 
galaxies and the intra-cluster gas. The considered interaction processes are galactic winds and ram-pressure 
stripping, which both transfer metal-enriched interstellar medium into the intra-cluster gas and hence increase its 
metallicity. We investigate the efficiency and time evolution of the interaction processes by simulated metallicity 
maps, which are directly comparable to those obtained from X-ray observations. We find that ram-pressure 
stripping is more efficient than quiet (i.e. non-starburst driven) galactic winds in the redshift interval between 1 
and 0. The expelled metals are not mixed immediately with the intra-cluster gas, but inhomogeneities are visible 
in the metallicity maps. Even stripes of higher metallicity that a single galaxy has left behind can be seen. The 
spatial distribution of the metals transported by ram-pressure stripping and by galactic winds are very different 
for massive clusters: the former process yields a centrally concentrated metal distribution while the latter results 
in an extended metal distribution. 
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1. Introduction 

The components of clusters of galaxies - galaxies and 
intra-cluster medium (ICM) - interact with each other 
. in many ways. There is more and more observational ev- 
idence that various types of processes are at work, which 
remove interstellar medium (ISM) from the galaxies. This 
metal enriched ISM mixes with the ICM, so that the cur- 
rently observed ICM is a mixture of both. Hence the met- 
als are a good tracer for the present and past interaction 
processes between galaxies and ICM. The new X-ray ob- 
servations by CHANDRA and XMM make it possible to 
measure not only profiles but 2D distributions of metals: 
metallicity maps (e.g. Schmidt et al. 2002; Furusho et al. 
2003; Sanders et al. 2004; Fukazawa et al. 2004; Hayakawa 
et al. 2004) showing that the metallicity has clearly a non- 
uniform, non-spherical distribution. 

Already many years ago supernova-driven galactic 
winds were suggested as a possible ISM transfer mech- 
anism (De Young 1978). Recent observations and simu- 
lations indicate that other processes can contribute con- 
siderably, as well (e.g. Hayakawa et al. 2004). One pro- 
cess gaining increasing attention is ram-pressure strip- 
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ping (Gunn & Gott 1972). In the Virgo cluster alone 
at least 7 ram-pressure affected spiral galaxies have been 
found (e.g. Cayatte et al. 1990; Vollmer et al. 2004). Also 
other processes like galaxy-galaxy interactions (Kapferer 
et al. 2005), jets from AGNs and intra-cluster supernovae 
(Domainko et al. 2004) can contribute to the metal en- 
richment. 

Currently, several approaches are being made to ex- 
plain the overall enrichment of the ICM by taking into 
account galactic winds. De Lucia et al. (2004) and 
Nagashima et al. (2004) use a combination of semi- 
analytic techniques and N-body simulations to calculate 
the overall ICM metallicity, but they do not predict the 
distribution of metals in a cluster. They find that mainly 
the massive galaxies contribute to the enrichment and that 
there is a mild metal evolution since z=l. Tornatore et al. 
(2004) do simulations with smoothed particle hydrody- 
namics that include detailed yields from type la and II 
supernovae, but do not distinguish between the different 
transport processes. They put most of their emphasis on 
the amount of iron produced and on metallicity profiles. 

Other groups have calculated the effect of ram- pressure 
stripping on single galaxies (e.g. Abadi et al. 1999; Quilis 
et al. 2000; Toniazzo & Schindler 2001; Roediger & Hensler 
2004), but not the effect on the ICM. We take a comple- 
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mcntary approach by investigating the transport of the 
metals from the galaxies into the intra-cluster medium 
and study in detail the efficiency of the various trans- 
port processes. In our simulations presented here, we have 
taken into account two different transport processes - ram- 
pressure stripping and galactic winds. 

2. Simulations 

2.1. Numerical method 

We use N-body and hydrodynamic techniques together 
with a semi-numerical phenomcnological galaxy forma- 
tion code and various prescriptions for interactions be- 
tween galaxies and the intra-cluster medium. The core of 
the programme is a hydrodynamic code with shock cap- 
turing scheme (PPM, Colella & Woodward 1984), mul- 
tiple grid refinement (Ruffcrt 1992) and radiative cool- 
ing. This hydrodynamic code is combined with an N- 
body tree code (Barnes & Hut 1986) - providing actual 
galaxy orbits - with constrained random fields as ini- 
tal conditions (Hoffman & Ribak 1991) as implemented 
by van de Weygaert & Bertschinger (1996). It also in- 
cludes an improved version of the galaxy formation code 
by van Kampcn ct al. (1999), which provides galaxy prop- 
erties. We use a ACDM cosmology = 0.7 ,fl m = 0.3, 
cr 8 = 0.93, and h = 0.7) with 64 3 particles. While the N- 
body code and the galaxy formation trace the whole evo- 
lution from the beginning the hydrodynamic covers the 
rcdshift interval between 1 and with time steps of about 
15 Myrs starting with the initial condition of hydrostatic 
equilibrium. The hydrodynamic simulation is calculated 
on 4 levels of nested grids centred at the cluster centre, 
the largest being (20Mpc) 3 and the smallest (2.5Mpc) 3 . 
The number of grid cells at each level is 128 3 . In addi- 
tion we use various prescriptions for the different inter- 
action processes which transport metal enriched gas from 
the galaxies into the ICM. These transport processes are 
calculated in the full simulation volume. 

2.2. Galactic winds 

The amount of matter ejected by galactic winds is calcu- 
lated with an approach initiated by Breitschwerdt et al. 
(1991). The algorithm includes thermal and cosmic ray 
driven winds in disk galaxies. The wind code requires sev- 
eral galaxy parameters like halo mass, disk mass, spin pa- 
rameter, scale length of the components, temperature dis- 
tribution of the ISM, magnetic field strength and gas den- 
sity distributions as well as stellar density distribution. We 
perform parameter studies and summarise their results in 
a look up table for about 1000 different disk galaxies. The 
database includes galaxies with halo masses of 1 x 10 10 
M Q - 1 x 10 12 M Q . In addition 10 different disk masses 
for a given halo mass (0.01 - 0.5 x halo mass) and 10 
different spin parameters for each halo-disk combination 
are calculated. As we want to investigate the contribution 
of galactic winds to the enrichment of the ICM of quiet, 



non-starburst disk galaxies (= steady-state winds) the pa- 
rameters of the ISM and the star formation rate (1 M Q /yr) 
are kept constant. Following this approach we find an av- 
erage mass loss rate of about 0.3 M /yr. In case of a mass 
loss due to a galactic winds the ISM is transferred with the 
corresponding metallicity and a temperature of 5 x 10 6 K 
into the hydrodynamic simulation at the position of the 
galaxy. 

2.3. Ram-pressure stripping 

Galaxies in galaxy clusters can suffer mass loss from inter- 
action with the ICM, e.g. ram pressure acts on the gas disk 
of galaxies moving through the ICM. We assume that gas 
is stripped off beyond the radius of the galaxy where the 
restoring gravitational force is equal to the force due to 
ram-pressure (stripping radius, Gunn & Gott 1972). We 
assume for disk galaxies an exponential stellar disk and 
gas disk. The validity of this approach has been shown by 
simulations and by comparisons of galaxies in the Virgo 
cluster affected by ram pressure (Abadi et al. 1999). For 
supersonic galaxies we follow the Rankine Hugoniot con- 
dition to derive the properties of the galaxies surrounding 
ICM. Mass loss of galaxies moving inclined through the 
ICM is scaled with the cosine of the inclination angle with 
respect to the mass loss derived for an uninclined galaxy. 
Galaxies which are stripped are assumed to have a trun- 
cated gas disk and they only lose more gas if increased 
ram pressure further reduces the stripping radius. 

3. Results 

We present simulations of clusters with different mass: a 
massive (1.3 x 10 15 M©, Cluster 1) and a less massive 
cluster (7.4 x 10 14 M , Cluster 2). Cluster 2 undergoes a 
major merger (mass ratio 1:3), whereas Cluster 1 only has 
small merger events. 

In Fig. 1 the X-ray emission weighted metallicity maps 
are presented for Clusters 1 and 2: only the enrichment 
that has taken place since redshift 1 has been taken into 
account. Cluster 1 is shown at redshift z = with two 
different enrichment mechanisms, ram-pressure stripping 
(Panel la) and galactic winds (Panel lb). The two mecha- 
nisms result in a very different spatial distribution. Ram- 
pressure stripping is very efficient in the centre, because 
there the ICM density is very high (~ 6 x 10~ 27 g/cm 3 
within a radius of 0.1 Mpc) and leads to a centrally con- 
centrated metallicity distribution. On the other hand the 
high central ICM density yields a high pressure which sup- 
presses galactic winds resulting in a very extended dis- 
tribution. Apart from the differences in the centre both 
processes yield pronounced inhomogeneities in the abun- 
dances. Single galaxies leave traces of metals in shapes 
of stripes behind them. Both the centrally concentrated 
metallicities and the general inhomogeneities are in good 
agreement with observations. The increase of metallicity 
(X-ray emission weighted) since redshift 1 averaged over 
the central cluster region within a radius of 400 kpc is 0.09 
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and 4 x 10 8 in solar units for the ram-pressure and the 
winds model, respectively. 

The situation looks very different in a less mas- 
sive cluster. In Cluster 2 the winds arc less suppressed. 
Therefore the winds contribute also to the central metal- 
licity. Inhomogeneities in the abundances do not disperse 
immediately (compare insets) even though the ICM in this 
cluster is stirred by a merger. The stripes are gradually 
spread out and at the end a roughly homogeneous region 
of high metallicity is present at the centre. A given inho- 
mogencity at t—1.7 Gyrs will be spread out on average 
over a 30 times larger volume at the end of the simulation 
(t=8 Gyrs), depending strongly on the local dynamics of 
the surrounding ICM, e.g. shocks. 

Non-starburst driven winds are less efficient in trans- 
porting gas than ram-pressure stripping in the time in- 
terval between redshift 1 and investigated here. During 
the non-merger phases of Cluster 2 ram-pressure stripping 
is about a factor of 5 more efficient than non-starburst 
driven winds. Figure 2 shows the evolution of the metal 
mass loss for ram-pressure stripping and quiet galactic 
winds in a merger cluster. The different efficiencies are 
clearly visible. The increase of metallicity (X-ray emission 
weighted) since redshift 1 averaged over the central clus- 
ter region within a radius of 400 kpc is 0.07 and 0.0035 
in solar units for the ram-pressure and the winds model, 
respectively. 

The mass loss due to ram-pressure stripping shows 
large fluctuations compared to galactic winds, because 
non-starburst driven winds have an almost constant mass 
loss over a long period of up to 5 Gyrs, whereas ram- 
pressure stripping removes huge amounts of gas from a 
single galaxy on much shorter time scales. During the two 
merger events (at 3 Gyrs and 5.5 Gyrs) the mass loss due 
to ram pressure is much higher because of the enhanced 
ICM pressure. In contrast non-starburst driven galactic 
winds are suppressed in galaxies located in the high ICM 
pressure region. Therefore the mass loss due to winds is 
slightly lowered during mergers. 



4. Summary and Conclusions 

Different enrichment mechanisms - ram-pressure stripping 
and galactic winds - yield different metallicity distribu- 
tions. 

— In massive clusters ram-pressure stripping provides a 
much more centrally concentrated distribution than 
galactic winds, because galactic winds can be sup- 
pressed in the cluster centre while ram-pressure strip- 
ping is most efficient there. 

— Both processes yield abundance distributions with lots 
of inhomogeneities in good agreement with recent X- 
ray observations. 

— The inhomogeneities are not dispersed immediately, 
but are gradually spread out. 
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Fig. 2. Evolution of the metal transfer from the ISM to 
the ICM in solar masses per year for model Cluster 2. 
The metal mass loss for ram-pressure stripping and for 
non-starburst driven galactic winds are shown. 

— The efficiency of gas transfer due to ram-pressure strip- 
ping is higher compared to non-starburst driven galac- 
tic winds in the redshift interval between 1 and 0. 

The comparison of these results with new detailed ob- 
served metallicity maps will yield insights into the vari- 
ous interaction processes between the ICM and the cluster 
galaxies. The knowledge on their efficiency and timcscales 
will provide information on the formation and evolution 
of galaxies and their hosting clusters. In the future we 
will extend the simulations to higher redshifts (z > 1) so 
that we will be able to compare the simulated metallicities 
directly with the observed metallicity maps. Further simu- 
lations including not only more interaction processes, but 
also distinguishing between different chemical elements 
are planned. 
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Fig. 1. X-ray weighted metal maps of two simulations at present epoch and insets at an epoch 6.7 Gyrs ago (a and 
b, Cluster 1: a massive galaxy cluster; c and d, Cluster 2: a merging cluster). Panels a and c show the metal maps for 
ram-pressure stripping and panels b and d those for galactic winds (non-starbursts). The quantities are given in solar 
abundances and the size of the maps is 5 Mpc on the side. Ram-pressure stripping always yields a centrally concentrated 
metallicity distribution. The distribution of metals transported by winds depends strongly on the cluster mass: in a 
massive cluster (b) the winds are suppressed in the centre resulting in a low abundance there. The comparison with 
the insets shows how the inhomogeneities spread/evolve. 



